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Abstract 

I The covariant particle-vibration coupling model within the time blocking approximation is em- 



ployed to supplement the Relativistic Random Phase Approximation (RRPA) with coupling to 
collective vibrations. The Bethe-Salpeter equation in the particle-hole channel with an energy de- 
pendent residual particle-hole (p-h) interaction is formulated and solved in the shell-model Dirac 
^ I basis as well as in the momentum space. The same set of the coupling constants generates the Dirac- 



Hartree single-particle spectrum, the static part of the residual p-h interaction and the particle- 



^ , phonon coupling amplitudes. This approach is applied to quantitative description of damping 

I phenomenon in even-even spherical nuclei with closed shells ^osp]-, ^^^^^ ^^^Sn. Since the phonon 

coupling enriches the RRPA spectrum with a multitude of phCgiphonon states a noticeable fragmen- 
tation of giant monopole and dipole resonances is obtained in the examined nuclei. The results are 
compared with experimental data and with results of the non-relativistic approach. 
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I. INTRODUCTION 



Recent development of experimental facilities with radioactive nuclear beams has stimulated 
enhanced efforts on the theoretical side to understand the dynamics of the nuclear many- 
body problem by microscopic methods. The most successful schemes based on the mean field 
concept use a phenomenological ansatz incorporating as many symmetries of the system as 
possible and adjust the parameters of functionals to ground state properties of characteristic 
nuclei all over the periodic table. Of particular interest are the models based on covariant 
density functionals because of their Lorentz invariance. A large variety of nuclear 

phenomena have been described over the years within this kind of models: the equation of 
state in symmetric nuclear matter, ground state properties of finite spherical and deformed 

u n n 

nuclei all over the periodic table p|] from light nuclei |^ to super-heavy elements [5|], from 
the neutron drip line, where halo phenomena are observed [^| to the proton drip line Ql 
with nuclei unstable against the emission of protons j^. In the small amplitude limit one 
obtains the relativistic Random Phase Approximation (RRPA) ^j. This method provides a 
natural framework to investigate collective and non-collective excitations of p/i-character. It 
is successful in particular for the understanding of the position of giant resonances and spin- 
or/and isospin-excitations as the Gamov Teller Resonance (GTR) or the Isobaric Analog 
Resonance (lAR). Recently it has been also used for a theoretical interpretation of low-lying 
dipole [l^ and quadrupole excitations. 

Of course the density functional theory based on the mean field framework cannot provide 
an exact treatment of the full nuclear dynamics. It is known to break down already in ideal 
shell-model nuclei such as ^osp]^ with closed protons and neutron shells. In self-consistent 
mean field calculations one finds usually a considerably reduced level density at the Fermi 
surface as compared with the experiment. The RRPA describes very well positions of giant 
resonances but underestimates their width considerably. To solve the level density problem 
the covariant theory of particle-vibration coupling has been developed and applied in the 
Ref. 

In the present work we formulate the covariant response theory employing the 



particle- vibration coupling model within the time-blocking approximation |l3, 



15 



Ifl to 



describe the spreading of multipole giant resonances in even-even spherical nuclei. 
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II. FORMALISM 



In the relativistic nuclear mean-field theory motion of a single nucleon is described by the 
Dirac equation with an effective mass m* and a generalized four-vector of the momentum 

(7'^P^-m*)|^) = 0. (1) 

These quantities are determined by the scalar and vector = (S°, S) parts of the mass 
operator (self-energy) S within mean-field approximation: 

m* = m + t^s , Pf, = p^, - tf, = [i^ - So, iV + fl^ , (2) 

Eg is generated by the scalar cr-meson field. When we go beyond the mean-field approxima- 
tion, we have to take into account that in the general case the full self-energy S is non-local 
in the space coordinates and also in time. This non-locality means that its Fourier transform 
has both momentum and energy dependence. Let us decompose the total self-energy matrix 
into two components, a static local and and an energy dependent non-local term: 

S(r, r'; e) = t{r)5{v - r') + S^(r, r'; e), (3) 

where index "e" indicates the energy dependence. Due to time-reversal symmetry and the 
absence of currents the space-like components of S vanish, therefore only scalar and the 
time-like components of the mean field are considered in the following. 
The quantity S(r) is assumed to be the RMF self-energy generated by the cr, uj and £-meson 
fields within the framework of the no-sea approximation (see, for instance, Ref jSj). To 
describe the non-local part S'^(r,r'; e) we apply the covariant version of the particle-phonon 
coupling model Due to the decomposition it is convenient to work in the shell- model 
Dirac basis (I'i/'fc)} which diagonalizes the energy-independent part of the Dirac equation: 

/i^'l^fc) =£fc|^fc), /i^ = ap + /5(m + !],)+ So, (4) 

where hP denotes the Dirac hamiltonian with the energy-independent mean field. In the case 
of spherical simmetry the spinor is characterized by the set of single-particle quantum 
numbers k = {{k),mk},{k) = {n^, j^, tt^, tfc} with the radial quantum number n^, angular 
momentum quantum numbers jki'^^k-, parity tt^ and isospin tk- In this basis the Dyson 
equation for the single-particle Green's function can be formulated as follows: 

Y,{{e - ek)6ki - T.li{e)]Gik'{e) = 5kk'. (5) 
I 
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Matrix elements of the energy- dependent part of the mass operator 

Sh(^) = / d'^rd'r' Vfc(r)S^(r, r'; e)Mr') (6) 



are expressed in terms of the particle-phonon couphng model 

g{o-„) g(o-„)* 

'fcra I In 

^ „ e-e„-Or„(fi9_^r7) 

through the phonon vertexes 7"? and their frequencies Q'^. They are determined by the 
following relation: 

^ki = y^^kv,ik'Pk'i'^ Vki',ik' = -r^- (^) 

Vku^ik' denotes the matrix element of the residual interaction which is a functional derivative 
of the relativistic mean field with respect to nuclear density p and ffl is the transition density. 
Here we use the linearized version of the model which assumes that p'^ is not influenced by the 
particle-phonon coupling and can be computed within the relativistic RPA. In the present 



work the residual interaction is generated by the relativistic NL3 Lagrangian [17[. In the 
Eq. (|7j) cr„ = +1 if r?, is an unoccupied state of p- or a-types and cr„ = — 1 for an occupied 
n state of /i-type, rj — > +0. a denotes states in the Dirac sea with negative energies which 
arise in the Lehmann expansion of the single-particle Green's function due to the no-sea 
approximation. 

Equation has been solved numerically in the Ref. in the shell-model of Dirac states. 
A noticeable increase of the single-particle level density near the Fermi surface relative to the 
pure RMF spectrum is obtained for ^°^Pb. This improves the agreement of the single-particle 
level scheme with experimental data considerably. For the four odd mass nuclei surrounding 
^°^Pb the distribution of the single-particle strength has been calculated and compared with 
experiment as well as with the results obtained within several non-relativistic approaches. 
The nuclear dynamics of an even-even nucleus in a weak external field is described by the 
linear response function which is a solution of the Bethe-Salpeter equation (BSE) in the 
particle-hole (p-h) channel. In the beginning it is convenient to consider this equation in the 
time representation. Let us include the time variable into the set of single-particle quantum 
numbers and use the following number indexation to simplify the expressions: 1 = 
In this notation the BSE for the response function R reads: 

R{IA, 23) = -(5(1, 3)G(4, 2) + - ^ ^(1, 5)G(6, 2)W^(58, 67)i?(74, 83), (9) 

^ 5678 




FIG. 1: Bethe-Salpeter equation for the p-h response function in the graphical representation. Solid 
lines with arrows denote one-body propagators through the particle, hole or antiparticle states, 
weavy lines denote phonon propagators, empty circles are the particle-phonon coupling amplitudes 
and the small black circle means the static part of the residual p-h interaction. 

where 

W{U, 23) = U{U, 23) + 3)G~\A, 2) + tG-\l, 3)S^(4, 2) - tJ:'{l, 3)S"(4, 2). (10) 

Here the summation over number indices implies also integration over respective time vari- 
ables. G is the mean field single-particle Green's function and U is irreducible in the p-h 
channel amplitude of the effective interaction including an induced interaction due to the 
phonon exchange. The graphical representation of the Eq. (jHl) is shown in Fig. H 
Using Fourier transformation of the Eq. one comes to an integral equation where both the 
solution and the kernel are singular with respect to energy variables. Another difficulty arises 
because the Eq. © contains integrations over all time points of the intermediate states. This 
means that many configurations which are actually more complex than lplh(8)phonon are 
contained in the exact response function. In the Ref. the special time-projection tech- 
nique was introduced to block the p-h propagation through these complicated intermediate 
states. It has been shown that for this type of response it is possible to reduce the integral 
equation to a relatively simple algebraic equation. Obviously, this method can be applied 
straightforwardly to our case. The full formalism can be found in the Ref. [l8|. Making 
use the above mentioned time projection one can transform the Eq. (JH)) to the following 
algebraic equation within the so-called time-blocking approximation: 

RkikiMkii^) = RkikiMksi'-^) 

~ ^ RkikeMksi^) [^ksksMkr ^ksksMkri'^)] R'krkiMkai'^) ^ (11) 
k^kekrks 
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where R is the mean-field p-h propagator, V is the residual interaction defined by the Eq. (jH}, 
$ is the particle-phonon coupling amplitude including phonon contribution both into the self- 
energy and into the induced interaction. It is supposed that the summation is carried out 
over the whole Dirac space. 

Thus, to describe the observed spectrum of the excited nucleus in the weak external field P 
within this formalism one needs to solve the Eq. (jllj) and to calculate the strength function: 

S{E) = ^ lim^lm P;^^^Rk,k,Mk.{E + iA)P,^u,. (12) 

The imaginary part A of the energy variable is introduced for convenience in order to obtain 
more smoothed envelope of the spectrum. This parameter has a meaning of an additional 
artificial width for each excitation. This width emulates effectively contributions from con- 
figurations which are not taken into account explicitely. 



III. RESULTS AND DISCUSSION 



The developed approach is applied to a quantitative description of isoscalar monopole and 
isovector dipole giant resonances in the even-even spherical nuclei ^°'^Pb and ^'^^Sn. Details 



of our calculation scheme are given in the Ref. [18|. First, the Dirac equation for single 
nucleons together with the Klein-Gordon equations for meson fields (RMF problem) are 
solved simultaneously to obtain the single-particle basis. Second, the RRPA equations 
are solved to determine for the above mentioned phonons. These two sets form the multi- 
tude of ph®phonon configurations which enter the particle-phonon coupling amplitude $. 
Third, an equation for density matrix variation (convolution of the Eq. (fTT|) with the ex- 
ternal field operator) is solved with this additional amplitude. It provides an enrichment of 
the calculated spectrum as compared to the pure RRPA. The equation for density matrix 
variation has been solved both in the momentum and in the Dirac spaces to ensure propri- 
ety of our calculational scheme and identical results have been obtained. The calculated 
strength functions for the isoscalar monopole resonance in ^osp]^ ^^^Sn computed within 
the RRPA and the RRPA extended by the particle- vibration coupling (RRPA-PC) are given 
in the Fig. El The fragmentation of the resonance caused by the particle-phonon coupling is 
clearly demonstrated although the spreading width of the monopole resonance is not large 
because of a strong cancellation between the self-energy diagrams and diagrams with the 
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FIG. 2: Isocalar monopole resonance in Pb and Sn obtained within two approaches: the 
RRPA (dashed hnes) and the RRPA with the particle-phonon couphng RRPA-PC (sohd hnes). 
Both calculations are based on the relativistic Hartree (RH) approach with the parameter set NL3. 

TABLE I: Lorentz fit parameters of isoscalar EO strength function in ^'^^Pb and ^^^Sn 
calculated within the RRPA and the RRPA extended by the particle-phonon coupling model 

(RRPA-PC) as compared to experimental data. 







<E> (MeV) 


r (MeV) 




RRPA 


14.13 


1.17 


208pb 


RRPA-PC 


14.02 


1.57 




Exp. [19] 


13.73(20) 


2.58(20) 




RRPA 


16.13 


1.96 


132Sn 


RRPA-PC 


16.07 


2.37 
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FIG. 3: Isovector El resonance in Pb: the results obtained within the non-relativistic approach 
(left panel) with Woods-Saxon (WS) single-particle input and Landau-Migdal (LM) forces and 
calculations performed within the covariant theory (right panel) based on the relativistic Hartree 
(RH) approach with the NL3 mean field parameter set. The RPA calculations are shown by the 
dashed curves, the RPA-PC calculations - by the thick solid curves. Experimental Lorentzian is 
given by the thin solid curves. 

phonon exchange (see Fig. QJ. The mean energies and widths of these resonances are pre- 
sented in the Table H] As experimental data we display the numbers adopted in the Ref. 
[lol for the calculation of the nuclear matter compressibility from the evaluation of a series 
of data obtained in different experiments for the isoscalar monopole resonance in ^''^Pb. 
The calculated photoabsorption cross sections for the isovector dipole resonance in ^osp^j 
and ^^^Sn are given in the Figs. 01 and ID respectively. The left panels show the results 
obtained within the non-relativistic semi-phenomenological approach developed in the Ref. 
jlG^ which includes particle-phonon coupling on the base of Woods-Saxon single-particle 
input and Landau-Migdal forces. The right panels show the results of the calculations 
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FIG. 4: The same as in Fig. El but for ^^'^Sn. 



TABLE II: Lorentz fit parameters of the El photoabsorption cross section in 208p]^ g^j^^j ^^^Sn 
calculated within the RRPA and the RRPA extended by the particle-phonon coupling model 

(RRPA-PC) as compared to experimental data. 







<E> (MeV) 


F (MeV) 


EWSR (%) 




RRPA 


13.1 


2.5 


120 


208pb 


RRPA-PC 


12.8 


3.8 


114 




Exp. [20] 


13.4 


4.1 






RRPA 


14.7 


3.0 


112 


132Sn 


RRPA-PC 


14.3 


3.8 


108 




Exp. [21j 


16.1(7) 


4.7(2.1) 
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within the relativistic approach developed in the present work. To make the comparison 
reasonable, calculations within the non-relativistic framework have been performed with 
box boundary conditions for the Schrodinger equation in the coordinate representation which 
ensures completeness of the single-particle basis. In both cases one can see the noticeable 
fragmentation of the resonances due to the particle-phonon coupling. Moreover, one can find 
more or less the same level of agreement with experimental data for these two calculations. In 
case of the isovector El resonance in ^^^Sn this is, however, not so clear because the integral 
characteristics of the resonance obtained in the experiment of tthe Ref. [2^ are given with 
relatively large discrepancies. But the difference which is of fundamental importance is that 
in the semi-phenomenological approach one usually fits the parameters on all three stages of 
the calculation: first, the Woods-Saxon well depth is varied to obtain single-particle levels 
equal to experimental values, second, one of the Landau-Migdal force parameters is adjusted 
to get phonon energies at the experimental positions (for each mode) and, third, another 
Landau-Migdal force parameter is varied to reproduce the centroid of the giant resonance. 
Although the varying of the parameters is performed in relatively narrow limits, in some cases 
it is necessary to obtain realistic results. In contrast, within the fully covariant microscopic 
approach developed in the present work no adjustment of parameters is made. 
The mean energies and widths of the isovector El resonance computed within the RRPA and 
the RRPA-PC are displayed in the Table IHl For the isovector El resonance as well as for 
the isoscalar EO resonance Lorentz fits have been made according to the method developed 



in the Ref. [2^ within the energy region between one and three neutron separation energies. 
For the El resonance in ^^^Sn the lower energy limit is slightly higher in order to separate 
the distinct group of pygmy states from the giant resonance. 



IV. SUMMARY 



The relativistic random phase approximation is extended by the particle-vibrational cou- 
pling model. The Bethe-Salpeter equation is formulated in the two-body basis of Dirac 
states. Amplitude of the effective interaction entering this equation contains the static part 
originating from the pure relativistic mean field as well as the energy-dependent part caused 
by the particle-vibrational coupling. The latter term has been considered within the time- 
blocking approximation in a fully relativistic way using the covariant form of the nucleon 
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mass operator. 

The developed approach is apphed to the computation of spectroscopic characteristics of 
nuclear excited states in a wide energy range up to 30 MeV for even-even spherical nuclei. 
An equation for the density matrix variation is solved in the Dirac space as well as in the 
momentum space. The particle-phonon coupling amplitudes of collective vibrational modes 
below the neutron separation energy have been computed within the self-consistent RRPA 
using the parameter set NL3 for the Lagrangian. The same force has been employed in giant 
resonance calculations for the static part of the effective p-h interaction. Therefore a fully 
consistent description of giant resonances is performed. 

Noticeable fragmentation of the isoscalar monopole and isovector dipole giant resonances 
in 208p]^ ^^"^Sn is obtained due to the particle-vibrational coupling. This leads to the 
appearance of a significant spreading width as compared to RRPA calculations. This is 
in agreement with experimental data as well as with the results obtained within the non- 
relativistic approaches [3,13]. 
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